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We consider a stochastic model for the directed motion of a solid object due to the rectification 
of asymmetric surface vibrations with Poissonian shot-noise statistics. The friction between the 
object and the surface is given by a piecewise-linear friction force. This models the combined 
effect of dynamic friction and singular dry friction. We derive an exact solution of the stationary 
Kolmogorov- Feller (KF) equation in the case of two-sided exponentially distributed amplitudes. The 
stationary density of the velocity exhibits singular features such as a discontinuity and a delta-peak 
singularity at zero velocity, and also contains contributions from non-integrable solutions of the KF 
equation. The mean velocity in our model generally varies non-monotonically as the strength of the 
dry friction is increased, indicating that transport improves for increased dissipation. 



I. INTRODUCTION 

Understanding noise-induced transport on micro- and 
nanoscales in the absence of an external force gradient 
has been an active research area in recent years. In 
particular, it has led to the notion of Brownian mo- 
tors, which are used to model, e.g. the propulsion mech- 
anisms of cellular motor proteins. Here, directed motion 
results quite generically from the combination of (i) non- 
equilibrium energy input, e.g. due to chemical or mechan- 
ical fluctuations, and (ii) a spatial or dynamical asym- 
metry. A large variety of Brownian motor models im- 
plementing these principles have been developed, where 
typically the fluctuations are rectified using a ratchet- 
shaped potential; see the review articles [U [2] and refer- 
ences therein. 

In related work, small objects such as hydrogels, 
droplets and micron sized particles have been moved over 
solid surfaces using asymmetric surface oscillations or 
noisy vibrations without a net force This kind of 

noise-induced transport relies on a ratchet effect in the 
velocity coordinate, which requires in addition to (i) and 
(ii) a nonlinear resistance of the moving object to the im- 
posed force. Such a nonlinearity can originate, e.g. from 
nonlinear friction or wetting hysteresis. In the case of 
a solid object moving over a solid surface, the nonlinear 
dry (Coulombic) friction between the two solids leads to 
hysteresis of the contact line, which effectively rectifies 
the fluctuations: dry friction makes the object stick to 
the surface so that a threshold force is needed to move 
it. Directed motion results when the forces do not over- 
come the threshold symmetrically, and this may happen 
even though the total force is zero on average. 

A simple model for noise- induced transport can thus be 
implemented by coupling asymmetric fluctuations to an 
object subject to nonlinear friction. In Ref. [TU], shot- 
noise fluctuations given by a two-state (dichotomous) 
Markov process were investigated with different types of 
nonlinear friction. In our recent letter |llj . we consid- 
ered singular dry friction in combination with one-sided 



Poissonian shot noise (PSN). PSN can be considered as a 
standard model for a mechanical random force and gener- 
alizes the usual Gaussian noise, to which it converges in a 
certain limit. Since PSN contains an intrinsic timescale 
(the mean waiting time between successive shots) and 
also allows for a variation of the shot amplitudes, it is a 
suitable model for a broad range of different fluctuation 
scenarios. 



The interplay of dry friction and PSN not only leads to 
a non-zero mean velocity of the object in the absence of 
a mean force input, but also to non-monotonic and sin- 
gular features in the transport properties, as discussed 
in Ref. [TT]. In the present work, we generalize these 
results and obtain in particular an exact solution of the 
Kolmogorov-Feller (KF) equation for the stationary ve- 
locity distribution for the case where the PSN is two- 
sided with exponentially distributed amplitudes. In this 
case, non-integrable solutions of the KF equation con- 
tribute to the stationary distribution in a certain param- 
eter regime. The mean velocity as a function of the dry 
friction strength varies non-monotonically as the dry fric- 
tion is increased, reaching a maximum at a particular dry 
friction strength. The transport in our model therefore 
generally improves for larger friction, over a considerable 
range; this is highly counterintuitive. 

The paper is organized as follows. In Sec. [IT] we intro- 
duce our model for the directed motion of an object mov- 
ing under the influence of dry friction and PSN. General 
results on the form of the stationary velocity distribution 
are obtained. In Sec. |III| we derive the explicit solution 
for PSN with a two-sided exponential amplitude distribu- 
tion. Singular features as well as contributions from non- 
integrable solutions of the KF equation are discussed. In 
Sec. [V] we investigate the non-monotonic behaviour of 
the transport properties of our model. Possible general- 
izations of our results and their relevance to experiments 
are discussed in Sec. EH 
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II. DIRECTED MOTION DUE TO DRY 
FRICTION AND PSN 

We consider the dynamics of a solid object moving over 
a solid surface subject to friction and noise In one 
dimension, to which we restrict ourselves for simplicity, 
and if we set the mass of the object to m = 1, its velocity 
satisfies the equation of motion 



v(t) = --v{t) - a(v)A + £(t) - a. 

T 



(1) 



The friction between the two solids is expressed phe- 
nomenologically by two friction terms: a dynamic fric- 
tion linear in v with strength 1/r, where r is the inertial 
relaxation time, and a dry friction — o~(v) A with strength 
A. The sign function a(v) is defined as 



<r{v) 



1, 



0, 



-1. 



v < 



v = 



v < 0. 



(2) 



The singularity in o~(v) at v = models the non-analytic 
behaviour observed in the dynamics of a solid object due 
to dry friction. In particular, the dry friction force in 
Eq. can lead to complicated nonlinear stick-slip dy- 
namics, where the object alternates in an oscillatory or 
chaotic way between sticking and sliding states [HI [T3] . 

The force a denotes a constant drift force. We take 
this as the mean of the noise setting a = (£(<)) so 
that overall there is no net force acting on the object. 
Without loss of generality we only consider a > 0, since 
the case of negative a follows by symmetry. Despite the 
absence of a net force bias, the object exhibits a non-zero 
mean velocity due to the rectification of the asymmetric 
noise. These transport properties will be discussed in 
more detail in Sec. IVl below. 

If we write F(v) for the sum of all the deterministic 
forces, Eq. can be written in the generic form 



v(t) = F(v)+m, 



(3) 



where F (v) denotes a piecewise-linear force with a dis- 
continuity at v = 0: 



F(v) = 



F_(v) 



v>0 



v < 0. 



The two branches of the force are given by 



1 



F±(v) = — («-u±), 

T 



(4) 



(5) 



where we define the two deterministic fixed points v± 
from the conditions F±(v±) = 0, 



v± = — (a ± A)r. 



(6) 



Below, in Sec. Ill A we consider a regularized version of 
F(v) where the discontinuity at v — is replaced by a 
linear interpolation over a small velocity range. 

We consider Poissonian shot noise (PSN) as a model 
for the noise f(t) in Eq. Q. PSN is a mechanical ran- 
dom force that can be considered a generalization of the 
usual Gaussian white noise. It can be represented by a 
sequence of delta shaped pulses with random amplitudes 
AM 



z(t)=J2A k S(t-t k ). 



(7) 



fc=i 



The waiting time between the successive pulses at times 
tfc is assumed to be exponentially distributed with pa- 
rameter A, where A is the rate at which pulses arrive. 
Then n, the number of pulses in time t, follows a Poisson 
distribution 



P(n) 



(\tr 



(8) 



with mean Xt. When a pulse occurs, its amplitude A k is 
sampled, independently for each pulse, from a distribu- 
tion p(A). The mean and the covariance of z{t) are then 
given by 



(z{t)) = \(A) i 



(9) 



(z(t)z(t'))-(z(t))(z(t')) = X(A 2 ) p 5(t-t'), (10) 

where (••■)„ indicates an average over the amplitude dis- 
tribution p(A). In the limit of vanishing pulse amplitudes 
A — > and a divergent pulse frequency A — > oo, with 
(A) p = and X(A 2 ) p = const. = 2D, PSN reduces to 

Gaussian white noise. 

In Eq. Q we choose as a superposition of two PSN 
processes 



£(t) = zx(t) + z 2 (t), 



(11) 



where z\(t) and 22 (i) are specified by the pulse rates Ai ; 2 
and amplitude distributions p\^{A). In order to have a 
noise with zero mean we set the constant drift force a to 



o = <£(*)> = A a (A) pi +X 2 (A) p2 , 
where the mean amplitudes are given by 



(A) 



A Pj {A)dA, 3 = 1,2. 



(12) 



(13) 



In the Gaussian limit of the PSN, Eq. converges to 
Brownian motion with dry friction, which has recently 
been investigated both theoretically and experimentally 
[T5H2"3"] . No directed motion results in this case unless 
an additional bias, such as a constant force, acts on the 
object. On the other hand, asymmetric PSN has been 
shown to induce directed transport for overdamped dy- 
namics in a symmetric potential |24| . The effect of a 
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two-state (dichotomous) shot noise was investigated in 
Ref. [10] for different types of nonlinear friction, but no 
singular features were reported. The particular case of 
one-sided PSN, discussed in Ref. [TT], is obtained from 
the present model by setting A 2 = 0. 

The time dependent velocity distribution p(v,t) obeys 
the Kolmogorov- Feller equation [521 HE] 



d_ 

Of 



p(v,t) 



d_ 
dv 

2 



F(v)p(v,t) 
5>i (<?(«- A, t)) 



Pi 



p(v,t)) ,(14) 



4=1 



where the average (p(v — A,t)) p indicates the convolu- 
tion 



(p(v-A,t)) p . = 



p[y- A,t)pj(A)dA. (15) 



We are interested in the properties of the stationary prob- 



ability distribution p{v). In the Gaussian limit, Eq. (14) 
becomes the Fokker-Planck equation with the piecewise- 
linear force F(v), 

^P(«»*) = -^(«)p(«.*) + -D^aP(».*)- ( 16 ) 

Under the conditions of stationarity and a zero prob- 
ability current, the stationary distribution p(v) is then 
obtained from 



F(v)p(v) 



D^piv) 



(17) 



which can be solved by integration. Eq. (17) indicates 
that p(v) is continuous even in the presence of the dis- 
continuous force: the discontinuity at v — on the left- 
hand side of Eq. ( 17) can be balanced by a discontinuous 



derivative on the right-hand side, so that the stationary 
distribution has a cusp singularity at the point where the 
force is discontinuous. 

In the following we focus on the non-Gaussian param- 
eter regime of the PSN. Here, the stationarity condition 
becomes, instead of Eq. (17), 



F(v)p(v) = / G(v — v')p(v')dv' , 



(18) 



where the diffusion operator is expressed in terms of a 
Green's function (see below). The integral equation (18) 
indicates that PSN induces a non-local diffusion of the 
degree of freedom, in contrast to the local normal diffu- 
sion that is generated by standard Gaussian white noise. 
Eq. ( 18 ) can be derived in a straightforward matter by 
performing a Fourier transformation (FT) of Eq. (14) 
with p = 0. In Fourier space, where we denote Fourier 
transforms by a ~, the Green's function is given by 



G(k) 



4=1 



Pj(k) - 1 



(19) 



This leads to the representation 



2 ~ v 

G(v) = Y, X 3 / Pj(v')dv'-{\i + \ 2 )6(v). (20) 

• 1 J —OO 



4=1 



This expression can also be obtained directly by inte- 
gration of Eq. (18); 9{v) = (1 + a(v))/2 indicates the 
Heaviside step function. Because of this contribution 
the Green's function always has a jump discontinuity 
at v = 0, even if the amplitude distributions pj(A) are 
smooth. 

Let us now investigate the properties of p(v) at the 
discontinuity of the force. We assume that the stationary 
distribution has the form 



p{v) 



ciP+(v), v>0 
c 2 p-(v), v < 0, 



(21) 



where the constants C\ and c 2 guarantee normalization. 
Substituting v = Q ± into Eq. (18) yields 



^+(0+) c 2P _(0-) 



F_(0" 



(22) 



because the convolution integral on the r.h.s. is contin- 
uous at v = 0. The last result together with Eq. Q 
implies that p(v) must have a discontinuity at v = 0. 
Once the functions p±(v) have been found, Eq. ( [22] ) to- 
gether with the normalization constraint J p(v)dv = 1 
gives two equations for the two constants c±, c 2 . 

However, the above reasoning can be applied only 
when F + (0 + ) and F_(0~) have the same sign. When 
they do not, Eq. (22) clearly cannot be satisfied. In this 



case, there has to be a finite change in probability mass 
at the crossover from v < to v > 0, which can only 
stem from a delta peak at v = 0. We therefore need a 
more general ansatz 



(23) 



cip+(v), 


v>0 


r Q s(v), 


v = 


c 2 p~{v), 


v < 



The r.h.s. of Eq. ( 18 ) now has a jump discontinuity 
— (Ai + X 2 )T 0(v) from the convolution of the step func- 
tion in G with the delta peak in p(v). Other than this the 
r.h.s. is still contin uou s, so comparing values at v = + 
and v = 0~ in Eq. ( 18 ) gives 



F + (0+)c lP+ (0+) = F_(0-)c 2 p_(0-) - AT , (24) 
where we have abbreviated A = Ai + X 2 . When Tq = 



this reduces to Eq. ( 22 1 as it should 



In order to determine the three unknowns Ci, c 2 , Tq we 



now require one further condition in addition to Eq. ( 24 ) 
and the overall normalization condition. This condition 
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can be obtained from the derivative of the stationarity 
condition Eq. (fl8~| 



F'{v)p(v) + F(v)p'(v) 

= E A W Pj(v-v')p(v')dv'-Xp(v), (25) 



which also follows directly from Eq. ( 14 ) . Again compar- 



ing terms across the discontinuity leads to 

a (n(o + )p+(o+) + F+(o+y + (o+) + Ap + (0+)) 
= c 2 (F'_(0-)p-(0-) +F_(0-y_(0-) + Ap_(0-)) 



-ro^(A ift (o + )-A ift (o-)) 



(26) 



and this provides the third condition on the unknown 
normalization constants. The main problem is thus to 
determine the solutions p+(v) and p_ (v) on either side of 
the discontinuity. This would in principle require solving 
the integral equation ([l8]) , which is a difficult task. In the 
following we pursue an alternative solution method: from 
Eq. ( 18 ) we derive a local ordinary differential equation 



that can be solved separately for v < and v > 0. This 
approach is feasible for amplitude distributions pj that 
are given as single exponentials, as we demonstrate in the 
next section. 



III. EXACT SOLUTION FOR EXPONENTIAL 
AMPLITUDE DISTRIBUTIONS 

Exact expressions for the stationary distributions 
p±(v) can be derived if the amplitude distributions are 
taken as single exponentials 



P2 



, (27) 



so that the noise in Eq. is two-sided PSN. Here the 
parameters A±,A 2 > relate to the mean amplitudes 
defined above as (A) = Ai and (A) p2 = —A 2 , and the 
mean of the two-sided PSN is 



a = XiAi - X 2 A 2 



Eq. |27|) leads to the FTs 



1 - iAik 



1 + iA 2 k' 



(28) 



(29) 



Using Eq. (19) the Greens function G{k) is in this case 
Ai^i \ 2 A 2 



G(k) = - 



1-iAxk l + iA 2 k' 



(30) 



Multiplying the FT of Eq. (18) by (l+iA 2 k)(l-iA 1 k) one 



thus obtains an ordinary differential equation of second 



order 



1+Ai 



dv 



l-A 2 ~)F(v)p(v) 
dv ' 



-ap^ + XA^—piv). 

dv 



(31) 



This ODE can be solved separately for v > and v < 0, 
and since F(v) is linear in each range, these solutions can 
be found using FTs. These solutions are identical to the 
ones obtained directly from a FT of Eq. ( 18 ) under the 
assumption of a purely linear force F(v) = F±(v). We 
note, however, that for general amplitude distributions 
solving Eq. ( 18 1 in this way would not yield the correct 
solution. Instead one would need to derive and solve the 
corresponding local differential equation. For the expo- 
nential distributions considered here, the solution is then 



P±( v ) = f(v-v±), 



(32) 



where the two fixed points v± are defined in Eq. (JsJ) . 
In terms of the shifted velocity v = v — v±, f(v) is the 
solution of Eq. (31 1 with the linear force F(v) = —v/t: 



exp^-ikv-ir J G(fc')dfc'|dfc, 

/OO 
exp {— ikv} (1 — iA\k)' 
-oo 

(1 + iA 2 k)- X2T dk. 



\ — AlT 



For the following it is convenient to define 



f*(y) = Q(±v)(±v) 



f±(v) = Q(±v)(±v) 



Air-lg-5/Ai 

^2T—l v/A 2 



(33) 

(34) 
(35) 



so that, if we discard irrelevant proportionality constants, 



m = f?(i)*f 2 (t). 



(36) 



The solutions p± (v) follow by mapping back to the orig- 
inal velocity via v = v — v±. They are thus given as 
the convolution (indicated by "*" above) of two Gamma 
distributions, shifted by v±. Both p±(v) are normaliz- 
able (integrable) on the respective domains (— oo,0) and 
(0, oo), which is a consequence of the use of FTs. 

However, Eq. p6|) does not provide the complete so- 



lution of Eq. (31). We can also find non-integrable (over 
the positive or negative real axis) solutions of Eq. (31) 
if we use the Laplace transformation (LT). These solu- 
tions can contribute to the stationary density because of 
the piecewise-linearity of F(v), as discussed in more de- 
tail below. Similar to the FT case, the LT solutions of 
Eq. (31 ) will be given as p±(v) = h(v—v±), with h(v) the 
solution of Eq. (31) for F(v) = —v/t as before. Solving 
Eq. (31 ) for v > with the condition h(0) = using LTs, 



one obtains in a straightforward way (denoting LTs with 
a hat) 



h(s) = (l + A lS )-^ r (A 2 s-l) 



-A 2 r 



(37) 
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which yields upon inversion 

h{v) = f+{i)*.f+{i). 
Likewise one finds for v < 



(38) 



(39) 



Setting v = v — v± leads to the LT solutions of Eq. (31 ) 
for the piecewise-smooth force Eq. Q. 

Here, it is crucial to distinguish between the cases A > 
a and A < a, which are equivalent to F + (Q + )F_(Q~) < 
and F + (0 + )F_((T) > 0, respectively (cf. Eq. Q). In the 
case A > a we have v + < and u_ > 0, so that the two 
LT solutions are given by 



Pi(v) — Q(v)h(v — v + ) 
P2(v) = ©(— v)h(v — V-). 



(40) 
(41) 



However, neither can contribute to the stationary density 
since they diverge for v — > ±oo. In this case the solution 
is thus given by Eq. (23) with the p±(v) of Eq. (32 1. 
It contains the three undetermined constants ci, C2, To, 
which can be found from the normalization constraint 



and the consequences Eqs. (24) and (26) of the station- 
arity condition. 

In the case A < a, on the other hand, we have both 
v+ < and V- < 0. This yields the solution p\{v), 



Eq. (40), and in addition 



Pz(v) = Q(—v)0(v-—v)h(v — v-) 
Pa(v) = Q(—v)Q(v — V-)h(v — v_). 



(42) 
(43) 



The solutions pi(v) and p?,(y) diverge for v — > ±oo and 
have to be discarded. Only p^{v) is normalizable on the 
domain (-co, 0) and can therefore contribute to the over- 
all solution for v < 0. In the case A < a one thus obtains 



the stationary distribution in the form (21) with 



P+( v ) = f( v - "+) 



(44) 



as before, and 



P-(v) = f(v — V-) + C30(u — v-)h(v — V-), (45) 

A separate normalization constant C3 appears here for 
the LT contribution. Because of this, an extra constraint 
is needed to fix all of C\, C2, C3, beyond normalization 
and Eq. (22). This is again provided by Eq. (26), but 
with ro = 0. 

We note finally that f(v) and h(v) only exhaust 
three out of the four possible convolutions of the form 



/1 (v) * fi (w)< The fourth combination fi(v) * f% («) 
does not appear above because it is undefined: the inte- 
gral defining the convolution always diverges. 



A. Linearly smoothed force 

In order to understand in more detail the emergence 
of the delta-peak for A > a for the discontinuous force 



F(v), it is helpful to consider a linearly smoothed repre- 
sentation of the sign-function 



a e {v) 




v > e 

—e<v<e 
v < — e 



(46) 



which obeys 



a{v) = lim a e (v) 



e-s-0 



(47) 



In the linearized region, the force F(v), Eq. 0), is thus 
A 



F e (v) 



— e < v < e, 



(48) 



for small e. Performing the above derivation of p±(v) 
using FTs with this linear force leads to the stationary 
distribution 



Pe(v) (X 



a 
A 



Ae/A-l 



(49) 



for — e < v < e and small e. One can see that p £ (v) 
has a divergence at v = — ae/A and this lies inside the 
smoothing region for A > a. At the same time, p e (v) 
has finite probability mass when integrated over [— e, e]. 
It therefore turns into a delta-peak in the limit e — > of 
a genuinely discontinuous force. 



IV. FEATURES OF THE STATIONARY 
DENSITY 

We now explore the properties of the stationary veloc- 
ity distribution p(v) in the various parameter regimes. 
For definiteness we will set the time scale r = 1 and 
the velocity scale A2 = 1. We also fix a > so that 
X1A1 > X2A2; the case a < follows by symmetry. 
One of our more surprising findings will be that for both 
A > a and A < a the stationary velocity distribution 
can exhibit characteristic double peaks. 

One important feature of p(v) is that the building 
blocks f(v) and h(v), defined in Eqs. (36), (38) and (39), 



depend only on the rate and amplitude parameters A 1; 
A2 and A\. The friction strength A only appears in the 
way these functions arc shifted, via v = v — v± , and in 
the conditions for the coefficients c\ , C2 and C3 or Tq with 
which they are weighted in p(v). It therefore makes sense 
to consider first these more elementary functions, and in 
particular f(v) from Eq. (36). This function always fea- 



tures in p(v) while h(v) appears only for A < a. 

If we go back to how f(v) was calculated above, we 
see that it represents the stationary velocity distribution 
of a particle subject to a linear restoring force — v/t and 
two-sided PSN whose mean has not been subtracted off. 
Three different plots of f(v) are shown in the inset of 
Fig. [l] to illustrate the range of possible behaviours. One 
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sees that f(v) is skewed, and this is generically the case 
unless Ai = A2 and A\ — A 2 . 

When both noise rates are larger than the inverse in- 
ertial relaxation time, Ai,A2 > 1/t = 1, f(v) is smooth 
with a single maximum. The position of this maximum, 
which we will call v m , shifts to the right as Ai increases 
and becomes broader at the same time. This is as ex- 
pected since A\ governs the strength of positive noise 
pulses. Given our convention a > 0, the lowest value of 
Ai we are allowing is fixed by A1A1 = X 2 A 2 , and gives 
the smallest v m for fixed Ai, A2 and A 2 . At that point the 
sign of v m is determined by the faster noise, i.e. v m > 
for Ai > A2. 

The above trend, of v m increasing with A\ and the 
peak broadening, holds also if only one of the rates is low 
(< 1) but the other high (> 1). In that case, however, 
the minimum is "pinned" so that it is always on the side 
of the faster noise. E.g. for Ai > 1, A2 < 1 one always 
has v m > 0. 

Finally if both rates are low, the maximum of f(v) is at 
v m = 0. When also the total rate is low, i.e. Ai + A2 < 1, 
the maximum turns into a power-law peak 

/(C) oc {v]^**- 1 (50) 

near v = 0. This singular contribution is generically 
present also for larger rates but then has less dramatic 
effects. For example for A x < 1, A 2 < 1 but total noise 
rate Ai + A2 > 1, the leading order decrease of f(v) from 
its maximum at v m = follows the power law ( |50[ ). In 
the case Ai = A2 = 1, where the exponenet of this power 
law becomes one, an explicit form for f(v) can be found 
as 

f(v) = A -i\ A -i (e(^)e- CMl + S(-v)e^) . (51) 

In agreement with the general discussion above this has 
a peak at v m — and departs from this with terms pro- 
portional to \v\, leading to a jump discontinuity in the 
derivative at v = Q. 

Regarding the function h(v), we see directly from its 
definition in Eqs. (38 1 and (39) that for small v it will 



grow as |w| Al+A2 1 , thus matching the singular term (1 50h 



in /(«). 

In our further discussion we now need to distinguish 
A > a and A < a. 



The case A > a 



For A > a, p(v) is given by Eq. (23): the density 
consists of two distinct parts for v > and v < 0, and 
exhibits a delta-peak at v = in addition to the discon- 
tinuity, as illustrated in Fig. [TJ An additional maximum 
arises when the maximum in p+(v) (or p-(v)) y which is 
located at v m + v + (or v m + u_, respectively), is "visible", 
i.e. is not cut-off by the discontinuity. The conditions for 
the appearance of these two possible maxima are thus 
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FIG. 1: (Color online) The stationary density p(v) for A > a, 



Eq. (23 1, for parameter values Ai = 1.5, A2 = 3, Ai = 2, 
A = 2. The partial densities p±(v) — f(v — v±) are indicated 
by dashed lines. The red crosses are results from a direct 
numerical simulation of Eq. |TJ using an increment method 
for PSN [27]. Inset: The function f(v), Eq. (36]), for Ax = 2 
and rates (Ai, A2) as indicated. For Ai + A 2 < 1, f(v) has a 
power law peak at the origin. 
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FIG. 2: Characteristic double peaks in p(v) for A > a. 
Parameter values: (a) Ai = 1.5, A2 = 7.5, A\ — 5, A = 1. (b) 
Ai = 20, A 2 = 2,Ai = 0.1, A = 0.1. 



and 



= v m - (A + a) > 



(53) 



A-a<0, 



(52) 



where v m denotes the position of the maximum in f(v) 
as before. Since V- > and v + < 0, these conditions 
cannot both be satisfied, so at most one maximum can 
be present. We can then summarize the conditions as 
A < I a — v m \. Because we are looking at the regime 
A > a, a maximum can then be visible in some range of 
A only if - in f(v) - it is at v m < or at v m > 2a. This 
means in particular that the power law peak in f(v) that 
appears when Ai, A2 < 1 is never visible as it is located 
at v m = 0. 

Fig. [2] shows two plots of p(y) for A > a, which re- 
veal the peak, the step discontinuity, and a maximum 
to both the left and the right of the peak. Whether the 
discontinuity at v = has a positive or negative sign 
is determined by the relative location of the maxima in 
p±(v), the skewness of f{v), and the normalization con- 
stants C\ , C2 and Tq in Eq. ( 23 1 . 
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FIG. 4: Characteristic peaks of the stationary density p(v) 
for A < a. Parameter values: (a) Ai = A2 = 0.8, A\ = 3, 
A = 0.6a = 0.96. (b) Ai = A 2 = 0.45, Ai = 3, A = 0.4a = 
0.36. 



FIG. 3: (Colors online) The stationary density p(v) for A < a, 
Eq. (451, for parameter values Ai = A2 = 1.5, A\ — 3, A = 
0.4a = 1.2. The three distinct contributions to the density 
are indicated by dashed lines. Inset: p(v) for larger friction 
parameter A = 0.8a = 2.4. 



We close this subsection with a discussion of the physi- 
cal intuition behind the results for the stationary velocity 
distribution p(v). The regime considered here is the one 
of strong friction. The particle then experiences a restor- 
ing force for small velocity v that changes sign at v = 0. 
In the absence of a noise pulse, the particle velocity will 
therefore return to zero in finite time. With PSN, there 
is always a nonzero probability that the time interval be- 
tween two noise pulses will be long enough for this to 
happen. Therefore one has a finite probability in steady 
state of finding the particle with zero velocity, as repre- 
sented by the <5-peak in p(v). The secondary maximum 
in p(v) is inherited from the maximum in f(v) and re- 
flects the effects of positive and negative force pulses not 
having the same rates or amplitudes. 



1A1+A2- 



B. 



The case A < a 



In this case the density is given by Eq. ( 45 1 : the delta- 



peak disappears, but p(v) contains an additional contri- 
bution from the LT solution as indicated in Fig. [3j This 
contribution converges to the delta-peak as A — > a as 
illustrated in the inset of Fig. [3] Similarly to the regime 
A > a, the density can contain an additional maximum. 
This requires that one of the conditions Eqs. ( 52|53 ) be 
satisfied, though that is not necessarily sufficient: the ex- 
tra contribution h(v — V—) may cover up the maximum, 
as happens in Fig. [3] 

We find that a maximum in p(y) for v < is generi- 
cally visible for small A/a, but as A becomes larger the 
LT contribution increasingly dominates and eventually 
makes p(v) monotonically increasing again. In a narrow 
parameter range both the maximum inherited from f(y) 
and the increase towards v = from h(v — v_) are visible, 
giving two distinct peaks inp(w), as shown in Fig. [4ji. For 



1 from (501, 



Ai + A 2 < 1 the left of these peaks becomes a power law 
divergence with the power law \ v- 
cf. Fig. [4|z». 

It is possible also to find parameter settings where a 
maximum in p(v) is visible on the v > side, but this 
occurs only in a rather narrow parameter range. This 
maximum also never becomes singular as it is inherited 
from a maximum at v m > (rather than v m — 0) in 

m 

For a physical interpretation of the results above, one 
needs to recall that in the absence of noise impulses the 
particle is subject to the force F±(v) = —(v — v±)/t 
depending on the sign of v. As both i> + and V- are 
negative for the weak friction case A < a considered 
here, the velocity of the particle will thus approach V- = 
— (a — A)r, exponentially in time. When noise rates are 
low there is then a high probability in steady state of 
finding the particle with velocity v close to u_ , and this 
causes the peak at or near v— in p(v). 



V. TRANSPORT PROPERTIES 

As discussed in the introduction, directed motion, i.e., 
a non-zero mean velocity of the object, results from the 
coupling of asymmetric fluctuations to a nonlinear re- 
sponse. If we focus on the stationary regime again, the 
mean velocity is given by 



vp(v)dv, 



(54) 



and is generally non-zero for Ai 7^ A2 or A\ 7^ A2. A sim- 
pler expression for (v) can be derived by taking directly 
the noise average in Eq. and considering the steady 
state, where (v) =0. This yields 



(v) = -At{(t(v)} 

= At ( / p(v)dv 



p(v)dv 



(55) 



Eq. ( 55 1 immediately tells us that a non-zero mean ve- 



locity results physically from the presence of [TT] 
(a) Inertia (non-zero r) 
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FIG. 5: (Colors online) The mean velocity given by Eqs. (551 
and (57 1. The cusp singularity at A = a is visible, (a) Plot 



of (v) for increasing values of A\ (> A2). Here Ai = A2 = 
1.5. (b) For larger Ai (> A2), the maximum in (v) shifts to 
A/a = 1. Here A 2 = 1, Ai = 2. 



(b) The influence of dry friction (non-zero A) 

(c) An asymmetric stationary velocity distribution 

A positive mean velocity is induced when the total prob- 
ability of observing a negative velocity is larger than the 
total probability of observing a positive one. This is 
somewhat counterintuitive, but is here a consequence of 
the particular form of p(v). Since the friction is sym- 
metric, the asymmetry in the stationary distribution can 
only be induced by an asymmetric noise £(i). The asym- 
metric fluctuations are thus rectified by the nonlinear dry 
friction and lead to directed motion. 

Note that Eq. (55) neglects any contribution from a 
delta-peak at v = and is thus valid only for A < a. In 
order to obtain the equivalent expression for A > a one 
needs the average of u(v) over a delta-peak, which is a 
priori undefined because a(v) is discontinuous at v = 0. 
To get a meaningful result one has to consider the lin- 
earized sign- function Eq. (46 1. Integrating the linearized 
part over [— e, e] using the distribution Eq. (49) leads to 
a contribution 



a e (v)p e (v)dv oc 



— e 

a 
A 



a 
A 



Ae/A-l 



v a 
~e + A 



dv 

Ae/A-l 



dv + 0(e). 
(56) 



As e — > the term of order e vanishes, while the integral 
in the first term becomes the area of the delta-peak. In- 
cluding the proper normalization, the contribution to the 
average of the sign-function from the delta peak in p(v) 
can be read off as — oTq/A. For A > a the expression 
corresponding to Eq. ( 55 ) is then 



= A7 



p(v)dv - 



-r 



p(v)dv .(57) 



We are interested in the behaviour of the mean ve- 
locity as the strength of the dry friction is increased. 
Since we assume a > 0, the mean velocity is also always 
> 0. Clearly, for A = we have (v) = and trans- 
port is purely diffusive in this case. On the other hand, 



for A a the dry friction dominates the dynamics, so 
that (v) 0. Therefore, the mean velocity behaves in 
a non-monotonic way and reaches a maximum for a cer- 
tain value of A/a, which is shown in Fig. [5] The exact 
location of the maximum is difficult to determine ana- 
lytically due to the complicated form of the stationary 
distribution, but we can get a good qualitative picture 
by numerical evaluation. 

The transport generally improves as a increases. For 
X2A2 = const, the maximum of the mean velocity be- 
comes larger upon increasing either Ai or A\ and its po- 
sition shifts towards A/a = 1. This value indicates the 
change between the two distinct shapes of the stationary 
velocity distribution, discussed in Sec. [IV] In the plots of 
the mean velocity in Fig.[5]this transition manifests itself 
in a distinct cusp singularity at A = a. Note that when 
varying Ai or A\, the value of a changes, so the curves 
in Fig. [5] are each for a different range of A. For large 
a optimal transport is achieved for the dry friction value 
A = a. Physically, this optimal value can be understood 
by noting that a acts as a constant negative force bias 
in the equation of motion ([I]). For A = a this negative 
drift is cancelled exactly by the dry friction, so that mo- 
tion induced by the fluctuations is predominantly in the 
forward direction. This transport behaviour is similar to 
the one observed for the special case of one-sided PSN, 
discussed in Ref. [TT], since for larger a the effect of the 
fluctuations in the backward direction becomes negligi- 
ble. In the one-sided case, the critical friction strength 
A = a indicates a transition between two distinct modes 
of motion: directed motion for < A < a and unidirec- 
tional motion for A > a. In the unidirectional regime, 
the dry friction completely filters out the negative drift 
and no negative velocity fluctuations appear. 

For smaller a the mean velocity exhibits a maximum 
for A > a, since in addition to the negative drift the 
fluctuations in the backward direction are more relevant. 
The maximum is then also less pronounced. For larger 
A > a the dry friction again starts to increasingly inhibit 
the transport. Overall, our model therefore exhibits the 
striking effect that increasing dissipation enhances the 
transport in a wide range of parameter values. 



VI. DISCUSSION 

We have investigated an exactly solvable nonlinear 
model for the directed motion of an object subject to dry 
friction and Poissonian shot-noise (PSN) with zero mean. 
The transport in this model behaves in a non-monotonic 
way as the dry friction strength is increased, exhibiting a 
maximum at a certain dry friction strength. The singular 
features of the stationary velocity distribution should also 
appear for more general piecewise-continuous forces sub- 
ject to PSN. This follows from our discussion in Sec. [TTJ, 
where the discontinuity and the peak at v = are de- 
rived for a general piecewise force. However, the partial 
solutions p±(v) would be more difficult to obtain when 
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the force branches F±(v) are each given by a nonlinear 
force. The same would be true if one wanted to go be- 
yond a one-dimensional description. It is also not clear 
how non-integrable solutions would contribute for more 
general nonlinearities, given that there might be multiple 
velocity fixed points. 

We were only able to solve the KF equation for noise 
amplitude distributions that are given as superpositions 
of single exponentials. In Fourier space, single exponen- 
tials lead to a Green's function in the form of a rational 
function. Inverting the stationarity condition Eq. ( 18 ) 



in Fourier space then leads to an ODE for the station- 
ary distribution. Due to the local nature of the ODE, 
the solutions p±(v) can be obtained by considering the 
dynamics separately on either side of the discontinuity. 
On each side the force is linear and thus p± (v) can be 
derived in a straightforward way. These solutions are 
in fact identical to the ones obtained directly from the 
integral equation (18) by restricting the force to F±(v) 
separately. However, this direct solution of Eq. ( 18 ) is 



not correct for arbitrary amplitude distributions due to 
the mixing of the dynamics across the discontinuity in 
Eq. |T8l). 



Our results might be relevant for coarse-grained mod- 
els of recently developed granular motors. These mo- 
tors are in principle realizations of the classical Feynman 
ratchet and pawl system in a granular gas. In the exper- 



imental work of Ref. [3H] , a chiral rotor of millimeter-size 
is immersed in a gas of rigid particles. By coating the 
sides of the rotor with different materials, a ratchet ef- 
fect can be realized, leading to directed rotations. On 
a coarse-grained level the influence of the granular gas 
at different particle densities, acting like random force 
impulses, could be modeled using PSN. This could take 
into account different spatial and temporal noise char- 
acteristics. Theoretical studies based on the Boltzmann 
equation including dry friction effects have already found 
singular features in the velocity distribution similar to the 
ones obtained here [25] , 

Our model could be implemented in a straightforward 
way in experimental setups similar to the ones used by 
Chaudhury et al pTM23| . In fact, vibrations with PSN 
statistics should in principle be easier to implement than 
Gaussian noise, since the fluctuations appear on a sepa- 
rate timescale. Dry friction has recently been shown to 
be relevant in particular to biological soft matter systems, 
which opens up the possibility of realizing the transport 
properties of our model in a nanoscale system [19 | [30Tt32] . 
The thermal environment would provide an additional 
noise source, but overall the qualitative transport fea- 
tures should remain unchanged since thermal noise does 
not exert any net force. The key issue would be to find a 
suitable realization of the PSN, e.g. through intermittent 
chemical reactions. 
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